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Abstract. We present a thorough investigation of stellar 
photometry based on HST imaging of crowded fields at 
85" and 10" from the centre of the high-surface bright- 
ness elliptical galaxy M32. The Principal Investigators 
of the present archive data have elsewhere presented an 
impressive colour- magnitude diagram of the field at 85". 
Based on the same data we enlarge on their photometric 
analysis and supplement with error estimators that more 
clearly show the implications of severe image crowding on 
the stellar photometry. We show that the faintest stars 
(/ <; 25.0, V <; 26.0) are found too bright by several tens 
of a magnitude. For the field at 10" we conclude that it 
is not possible to obtain reliable stellar photometry, stan- 
dard deviations being larger than 0.4 mag. Artificial-star 
experiments show that only very few of the brightest stars 
of the luminosity function can be expected to represent 
single objects, the majority being either spurious or not 
as bright as measured. The paper as such introduces and 
demonstrates basic guide lines which may be used when 
dealing with stellar photometry of severely crowded fields. 

Key words: methods: data analysis - galaxies: elliptical 
- galaxies: individual: M32 - galaxies: stellar content 



1. Introduction 

In order to increase our knowledge of extra-galactic stel- 
lar systems like elliptical galaxies, we observe the aver- 
age galaxy brightness, which gives rise to the traditional 
broad-band colours and spectral- line indices. In addition, 
we observe the variance of the galaxy brightness, i.e., the 
surface-brightness fluctuations (Tonry & Schneider 1988). 
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Of course, the ultimative goal is to obtain stellar luminos- 
ity functions and colour-magnitude diagrams, i.e., accu- 
rate photometry of individual stars. 

Until recently, colour-magnitude diagrams have been 
obtained solely from ground-based imaging of stars in our 
local neighbourhood, i.e., field stars and stars in open- 
and globular clusters of our Galaxy, in addition to ap- 
proaches to Local Group dwarf spheroidals and the Magel- 
lanic Clouds (low surface-brightness stellar systems) from 
ground (e.g. Freedmann 1992) and space (e.g. Mighell & 
Rich 1996). Then Grillmair et al. (1996) presented an im- 
pressive colour-magnitude (CM) diagram of M32, our clos- 
ests elliptical and companion to the Andromeda galaxy, 
based on Hubble Space Telescope imaging (archive data, 
ID 5233). Using the same data we have been able to re- 
produce the CM diagram of M32 which in itself, however, 
was not the main goal of our work. We have carried out an 
investigation of the reliability of the CM diagram, that is, 
an estimation of the photometric quality of the detected 
stars. The image crowding in dense fields like that of the 
high surface-brightness elliptical M32 may cause spurious 
detections and severely reduce the reliability of the pho- 
tometry. In the present paper we discuss specific aspects of 
the photometric quality of extra-galactic stellar systems, 
aspects that supplement the discussion by Grillmair et al. 

In our latest paper (Sodemann & Thomsen 1996) we 
presented ground-based imaging of M32 obtained with the 
Nordic Optical Telescope with a seeing of FWHM ~ l'.'O 
using 0"175 pixels. We identified a systematic variation in 
the I and S-band Surface-Brightness Fluctuations (SBF) 
of 0.2 — 0.3 mag in the radial range 10" ^ R ^ 50". In our 
search for the stars responsible for this variation we made 
a comparison of simulated stellar fields with the ground- 
based SBF observations. From this it became clear that it 
was not possible to supplement the measurements of the 
SBF with stellar photometry of the most luminous stars, 
due to excessive crowding. 

However, at least in some cases the superb resolution 
of the Hubble Space Telescope's (HST) Wide Field Plan- 
etary Camera 2 - the Planetary Camera (PC) with see- 
ing of FWHM O'.'l using 0'.'05 pixels - makes it possible 
to resolve more than just the very brightest stars. Grill- 
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mair ct al. (1996) presented a CM diagram based on HST- 
imaging of a region 85" from the centre of M32. From the 
surface brightness /U of M32 we can compare this off-centre 
field with a central field at 10" from the centre. Assuming 
that the same type of stars are responsible for the surface 
brightness at the two distances, the number of stars per 
area is ~ 40 times higher in the central field {ji 1 „ ~ 19.8 

and n 1 „ ~ 15.8). This merely indicates differing prob- 
lems introduced by image crowding, problems that may 
turn up when dealing with CM diagrams of fields at var- 
ious distances from the galaxy centre. We shall return to 
this later. 

When discussing luminosity functions and CM dia- 
grams we are interested in knowing to what magnitude 
limit we have actually been able to observe stars, a limit 
set mainly by a compromise between the image crowding 
and resolution, and the amount of exposure time. This 
magnitude limit is usually measured by addition and re- 
trieval of a small number of faint artificial stars, that is, the 
test of (in) completeness. However, whereas the traditional 
method of adding artificial stars estimates the detection 
probability it does not measure the photometric quality 
of the detected stars or the effects of spurious detections 
caused by severe image crowding. That is, the estimate 
of the limiting magnitude is based on information about 
the location (the coordinates) of the added/recovered star 
alone and not on the magnitude of the recovered star. For 
the data presented here, one important result of the image 
crowding is the following: An artificial star of known mag- 
nitude added to the programme frame may be recovered 
with a magnitude up to l m brighter than its original mag- 
nitude, but it turns out that it is not equally likely that 
this star will be recovered with a magnitude l m fainter 
than its original magnitude. (In Sec. 3 we show that the 
asymmetry in the distribution of recovered stars cannot 
be explained as a result of plotting the data as a function 
of magnitude). Thus, the typical measure a, the standard 
deviation, of the accidental error at a given magnitude is 
not an adequate description for this asymmetric distribu- 
tion of recovered stars and has to be supplemented with 
additional error estimates. 

Drukier et al. (1988) raised the problem of 'bin jump- 
ing' and pointed out a way to account for this when cor- 
recting a measured luminosity function for incomplete- 
ness. LePoy et al. (1993) also discuss the effects of severe 
image crowding, which may produce an artificial enhance- 
ment of the bright end of the measured luminosity func- 
tion if not dealt with in an appropriate way (see also Ren- 
zini 1992). 

We shall concentrate on the implications of crowding 
on stellar photometry of HST imaging of extra-galactic 
systems like the compact elliptical galaxy M32. We begin 
with a description of the data and image handling (Sec. 2). 
After that follows a discussion on the artificial-star exper- 
iments in connection with the off-centre field at 85" from 
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Fig. 1. Luminosity functions for objects identified in I and V, 
respectively, in the HST PC-frame at 85" from the centre of 
M32. A total of 19461 stars were found to match in I and V. 
No correction for M31 has been made, see text. The slope of 
the I-band luminosity function is found to be 0.42. (This value 
is not used elsewhere in this paper) . 

the centre of M32 (Sec. 3) and the dense central field at 
10" (Sec. 4). We summarise the results in Sec. 5. 

2. Data and image handling 

The observational basis for this paper is HST imaging of 
an off-centre field at 85" from the centre of M32 taken with 
the WFPC2 on October 22, 1994 (standard-calibration 
observations, archive data, ID 5233). See Grillmair et al. 
Fig. 1&2 for an outline of the M32 PC field on the sky. 
The 2x4 500 sec frames in F555W (V) and F814W (I), re- 
spectively, were co-added through the cosmic-ray remov- 
ing software of Groth & Shaya (1992) resulting in two 
programme frames. An I-band programme frame of the 
field at 10" was produced in a similar way from 3 expo- 
sures of 400 sec each (archive data, ID 5464). Correction 
for vignetting was carried out according to Holtzman et al. 
(1995a), and standard magnitudes and colour were derived 
using the calibrations in Holtzman et al. (1995b). We have 
adopted Aj = 0.14 for the galactic absorption (Burstein 
& Heiles 1984, Cardclli et al. 1989): The luminosity func- 
tions and CM diagrams have been corrected for this. 

The programme frames did not contain sufficiently 
bright, isolated stars to define a proper point-spread func- 
tion (PSF). Therefore, model PSFs generated by the 
program Tiny Tim (Krist 1993) defined the artificial 
stars added during the completeness tests and the mas- 
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Fig. 2. Colour-magnitude diagram of the off-centre field in 
M32. This figure was produced by matching the objects in 
Fig. 1. The 2x12 star signs show the magnitude and colour of 
the stars in our artificial-star experiments. 

tcr PSF of the psf-fitting software by Stetson's (1987) 
DAOPHOT II package. 

The completeness tests were carried out by adding 
stars distributed according to the observed CM diagram. 
The 2x12 star signs in Fig. 2 for I and V, respectively, out- 
line the range of magnitudes and colours. A total of 7110 
artificial stars were added to each of the two programme 
frames. When carrying out an addition of stars, only very 
few (in principle one, or several laid in a grid) ought to be 
added at a time so as to avoid confusion between added 
stars and an increase of the image crowding. As an alter- 
native, we split the addition of the artificial stars into two 
and those stars which happened to be added near another 
artificial star were removed. The stars were added using 
our HST-Sim routines, which scale and add the required 
number of PSFs (the routines are implemented using In- 
teractive Data Language, IDL). The routines take into 
account the geometric distortion inherent in the WFPC2 
optical system, and shifting and rebinning of the subsam- 



plcd PSFs provide subpixel positioning. Finally, Poisso- 
nian photon noise and Gaussian read noise are added. 

The photometry of objects in the programme frames 
was obtained through a procedure that included several 
passes of object identification and measurement of magni- 
tudes, using DAOPHOT II. First we worked through the 
identification of objects in the I-frame. In order to iden- 
tify as many objects as possible and to ease the identi- 
fication of faint objects we carried out ~ 6 iterations, 
successively reaching fainter magnitudes. For every two 
iterations we concatenated the lists of hitherto revealed 
stars, subtracted those from the original frame, and con- 
tinued from there. The iterations ended when we reached 
a threshold of 4 standard deviations of the background 
level, as calculated by DAOPHOT - at that point only 
a few hundred stars were identified and successfully sub- 
tracted by ALLSTAR and an additional iteration did not 
provide more objects. The final ALLSTAR list of all stars 
in I was then applied for star-subtraction in the V frame. 
After having calculated flux-averaged positions based on 
the I and V list of objects we again ran ALLSTAR on 
both frames, now with fixed coordinates. In this way a 
total of 19461 stars were found to match in V and I. The 
photometry of the star-added frames was carried out in a 
way identical to that of the programme frames. 

Note, Grillmair et al. applied Stetsons's (1994) pho- 
tometry package ALLFRAME. They found that the 
global characteristics of the CM diagram produced us- 
ing DAOPHOT II/ ALLSTAR were essentially identical to 
those found using ALLFRAME. 

We have made no attempts to correct for the M31 
background stars as none of the conclusions of this paper 
depend on knowing this accurately. According to Grillmair 
et al. about 18% of the ~ 20000 stars identified in the 
programme frames originate in M31. 

3. Results from the off-centre field 

The luminosity functions and CM diagrams of the off- 
centre field of M32 are seen in Fig. 1 and 2. With the 
artificial-star experiments as our basis we will now discuss 
the uncertainty associated with the identified stars. 

In Fig. 3 and 4 the magnitude of each group of added 
stars is indicated by a vertical line. The magnitudes of 
the recovered stars are binned into histograms of 0™1 
bins, one histogram for each magnitude group of added 
stars. Each of the histograms has been normalised so that 
they all cover the same number of stars (implying that the 
number- axis is arbitrary). This illustrates that the bright- 
est stars are easily identified, whereas fainter stars are 
'spread out' from their original magnitude (the vertical 
line) to surrounding magnitudes. Firstly, note that for the 
faintest bins the probability of a star being found in a bin 
other than the one to which it was added is nonvanishing. 
If we want to correct the luminosity function for incom- 
pleteness, this 'bin jumping' has to be taken into account 
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(Drukier et al. 1988). Secondly, the 'bin jumping' is clearly 
asymmetric with an extended tail of bright measurements. 
As the measurement errors are clearly not Gaussian we 
shall replace the usual, but now inadequate, standard de- 
viation a by the mean absolute deviation (MAD). (The 
MAD is related to the median as the standard deviation 
is related to the mean, but the MAD does not require 
that the distribution of the errors is Gaussian and it is 
a more robust estimator than is the mean, Press et al. 
1992). In addition, we supplement with the lower and up- 
per quartile, and several percentiles, of the magnitude of 
the recovered star m roc . These numbers are indicated in 
Fig. 3 and 4 and illustrate that the faintest stars in the 
CM diagram are very likely identified too bright by several 
tens of a magnitude as compared to their true, unaffected 
magnitude. 

Note, the straight lines represented by diamond sym- 
bols in Fig. 3, 4, and 9 indicate the magnitude of the 
added stars m a dd and the expected magnitude of the re- 
covered stars if they are not affected by 'binjumping', i.e., 
m rcc — m a( jd- In order to claim that we are able to carry 
out reliable photometry, m a dd must be a unique function 
of TO roc , and we must be certain that we are measuring the 
added star and not a statistical lump of unresolved stars. 

The distribution of the colour of recovered stars is, 
however, fairly symmetric and much less wide than that of 
the two passbands individually (Fig. 5). That is, the added 
stars which undergo 'bin jumping' tend to keep their orig- 
inal colour, mainly because they are most likely identified 
on top of a star with a colour similar to their own - the 
majority of stars in the CM diagram are found along the 
giant branch with (V — I) ~ 1.0. We may therefore ex- 
pect that the colours of the stars in the CM diagram are 
correct. 

Fig. 6 illustrates that the asymmetric distribution in 
the histogram of recovered stars is not a result of plot- 
ting the number of stars as a function of magnitude. We 
generated a Gaussian distribution of fluxes according to 
dN = Noe-V-hf/^df, where N = 10 6 , / = 1.0, and 
a = 0.4. The corresponding magnitudes m = —2.5 log/ 
were binned into 0™1 bins. The increase in the number of 
stars fainter than / (m > 0) exceeds a small increase just 
below m — and is clearly opposite the asymmetry seen 
in Fig. 3 and 4. Note that Seeker & Harris (1993) adopt a 
Gaussian function to represent the measurement error as 
a function of magnitude. 

The limiting magnitudes at half probability are V\ /2 = 
26.2 and I 1/2 = 25.1 for the M32 off-centre field. The limit- 
ing magnitudes are based on all the stars plotted in Fig. 3 
and 4, respectively, i.e., regardless of the associated error 
estimate. Note, the standard deviation a, as calculated by 
ALLSTAR, amounts to only m 12 and 0™16 at the limiting 
magnitudes for I and V, respectively, whereas the MAD is 
higher by a factor of 3.8 and 2.8 as compared to the stan- 
dard deviation (for a Gaussian distribution MAD ~ 0.8a). 
This merely indicates that in the present investigation the 
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Fig. 3. Results of our I-band artificial-star experiments. For 
each magnitude group of added stars (vertical line) the his- 
togram of recovered stars is shown. Note, the number-axis is 
arbitrary (Sec. 3). To clarify the implications of image crowd- 
ing we plot several error indicators, in addition to the mean 
absolute deviation (MAD) , of the magnitude of the recovered 
star m rcc . This demonstrates that the distribution of recov- 
ered stars is skewed towards brighter magnitudes, that is, the 
faintest objects are found too bright. This 'bin jumping' was 
discussed by Drukier et al. (1988). 



standard deviation, as given by ALLSTAR, may be con- 
sidered too optimistic. 



4. The M32 centre 

Based on the investigation of the off-centre field it is pos- 
sible to predict the number of identified stars in a central 
field at about 10", where the density is ~ 40 times higher 
than in the off-centre field (Sec. 1). The I-band luminosity 
function of the off-centre field starts to turn over at a num- 
ber density of ~ 1000 per (f^l in the PC frame (Fig. 2). 
In a field with density 40 times higher than this, the turn- 
over would correspond to a magnitude of ~ 20™5, or at 
the very brightest end of the luminosity function where the 
number density drops to zero. An interesting question is 
whether it is possible to find stars of 20™ 4 in the central re- 
gion of M32. (This magnitude corresponds to Mj ~ —4.0, 
the expected tip of the red-giant branch for metal-poor 
stars, e.g. Lee et al. 1993, and could provide an estimate 
of the distance to the galaxy). To address this question 
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Fig. 4. Same as Fig. 3 but now for the V-band. 



we have investigated a central field in M32 (archive data, 
ID 5464), the analysis being virtually identical to that of 
the off-centre field. (The addition of artificial stars was 
split up so that no more than one magnitude group at a 
time, corresponding to ~ 700 stars, was added to the pro- 
gramme frame). The luminosity function of the central 
field is seen in Fig. 7. 

Apparently we do see stars of magnitude 20^4 (Fig. 7). 
However, when we discuss the central field, we also have 
to consider the confusion that spurious detections due to 
SBF may cause, i.e., the confusion between single objects 
and detection of statistical lumps in the background of 
unresolved stars (they both have the shape of the PSF). 
This problem is illustrated in Fig. 8 where we see that 
bright objects around, e.g., 20^4 are indeed detected and 
located in the recovery tests with artificial stars of magni- 
tude 21™0 and fainter. The completeness tests show that 
if stars of magnitude 20^4 are present in the centre of 
M32, then we will be able to identify them (the data in 
Fig. 8 imply a limiting magnitude at half probability of 
Ji/2 = 21.6 for the central field of M32). What we do not 
know is how large a fraction of the supposedly identified 
bright stars in Fig. 7 are merely compounds of fainter, un- 
resolved stars. The limiting magnitude itself gives us no 
information about that. 

As will be evident from the following, we can only ex- 
pect very few of the brightest stars of the observed lu- 
minosity function to be single objects. Fig. 9 illustrates 



2.51.51.51.01.01.01.01.01.0 V- 



0.10 0.05 0.06 0.08 0.10 0.14 0.17 Q.23 0.22 MAD in mags 



-1 



20 22 



24 



26 



28 



Fig. 5. Colour deviation of recovered stars for the 9 artificial 
I-band magnitudes (error indicators are as in Fig. 3 and 4). 
This shows that stars are very likely to be found with the 
correct colour. The explanation is that the majority of stars 
in the CM diagram (Fig. 2) have more or less identical colours 
(V-7) ~ 1.0. 
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Fig. 6. Asymmetry that arises due to the effect of plotting a 
histogram of stars as a function of magnitude m = —2.5 log/, 
where the distribution in flux / is Gaussian (arbitrary num- 
ber-axis). The asymmetry is opposite that of the histograms 
in Fig. 3 and 4. 



that it is not possible to distinguish between the distri- 
butions of stars added at ~ 20™5 and fainter. That is, 
when we identify a star with magnitude ~ 20™5 or fainter, 
then we can by no means tell the true magnitude of that 
star. The completeness tests imply a limiting magnitude 



of h 



/2 



21.6 as mentioned above, but this limit only 



indicates the detection probability. We are inclined to dis- 
regard the stars in Fig. 7 that are fainter than ~ 20™5 (i.e., 
even stars brighter than Ii/ 2 ), since we have no reason to 
believe that they are single point sources. 

This is further emphasised by a simple simulation of 
the stellar field at 10": As the stellar content we adopt only 
one type of stars, namely stars of magnitude fhj = 22.78 
(Sodemann & Thomsen 1996). Our HST-Sim routines add 
the required number of such 'fluctuation stars' (PSFs) of 
magnitude mj to reach the average surface brightness of 
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Fig. 7. Luminosity function for objects identified in I in the 
HST PC-frame centered on M32. A total of 1606 objects were 
found in an annulus of area 80" and radius ~ 10". The median 
indicated at the top is 21™3. Only a few of the very brightest 
stars can be considered to be single objects. 



fj, „ ~ 15.8 as observed for the central frame, and pho- 
tometry is carried out by DAOPHOT II. The luminosity 
function of this simulated stellar field is seen in Fig. 10 and 
should be compared with the luminosity function of the 
observed central field (Fig. 7). The similarities are striking 
and tell us that the objects in Fig. 7 may all, except for 
a very few of the brightest, be described as mere fluctu- 
ations and not as single point sources. Based on Fig. 10 
we may explain the observed 'skewness' in Fig. 9 as fol- 
lows. Those of the faint added stars which were located on 
top of a bright fluctuation have inevitably been measured 
too bright and therefore end up in the left part of the 
'magnitude bands' of recovered stars (a magnitude band 
being defined by the 93.75 and 6.25 percentiles), whereas 
the stars in the right part of the 'magnitude bands' were 
located on top of a less bright fluctuation. 

We do not correct the observed luminosity functions 
for incompleteness and binjumping for the following rea- 
sons. Concentrating on the /-band, for the field at 85" 
the histograms of recovered stars are distinguishable for 
magnitudes that correspond to the giant branch and down 
to / ~ 24.0 where incompleteness starts to set in (Fig. 1 
and 3). In this paper we do not require the incompleteness- 
corrected luminosity function, and correction for binjump- 
ing would only affect the magnitudes around the limiting 
magnitude. For the field at 10", the luminosity function is 
comparable to any of the distributions of added test stars 
fainter than ~ 20™ 5, and an attempt to deconvolve the ef- 
fects of 'binjumping' would thus be an ill-posed problem. 



5. Conclusion 

In this paper we have discussed the effects that severe im- 
age crowding has on stellar photometry. The investigation 
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Fig. 8. Results of our I-band artificial-star experiments, now 
for the central region of M32, approximately 10" from the cen- 
tre. Again, the number-axis is arbitrary (Sec. 3). 



covered two fields of the high-surface brightness elliptical 
galaxy M32 observed with the HST's PC (archive data, 
ID 5233 and 5464). 

We have carried out a large number of artificial-star ex- 
periments, that is, addition and retrieval of artificial stars. 
Our experiments clearly show the presence of 'binjump- 
ing': The faint stars added to the programme frames are 
recovered too bright; the fainter the star and the denser 
the field, the bigger is the effect of 'binjumping'. The tra- 
ditional indicator of (in) completeness of the photometry, 
777,1/2, merely provides a detection probability but not the 
quality of the photometry. 

For the less dense field at 85" from the centre of M32, 
the effects of image crowding show up around the limiting 
magnitude ii/ 2 (25™1). Stars fainter than that are found 
too bright by several tens of a magnitude as compared 
to their true magnitude, whereas their colour is mainly 
unaffected. For the far more dense field at 10" from the 
centre, 'binjumping' affects stars that are even brighter 
than 1-1/2 (21™6). As our Fig. 9 demonstrates, only the 
very brightest stars with I 20'. n 5 may be considered to 
represent single objects. Thus, we do measure single stars 
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Fig. 9. Same data as Fig. 8 (besides three extra sets of artificial 
stars) but plotted in a slightly different way (error indicators 
are as in Fig. 3 and 4). As emphasized in the text, m a( jd must 
be a unique function of m roc , and we must be certain that 
we are measuring the added star and not a statistical lump 
of unresolved stars, before we can claim that we are able to 
carry out reliable photometry. This is clearly not the case for 
magnitudes fainter than ~ 20^5. From that magnitude and 
onwards it is no longer possible to distinguish between the 
distributions of recovered stars. Therefore, only a few of the 
very brightest stars with / ^ 20.5 in Fig. 7 may be expected 
to be single objects. 



of, e.g., 20v n 4 (corresponding to Mj ~ —4.0, the expected 
tip of the red-giant branch for metal-poor stars) in the 
off-centre field at 85", but based on the present data it 
is not possible to claim that for the central field at 10", 
where the density is 40 times higher. 

Finally, we would like to emphasize that the data 
analysis, which has been introduced and described in the 
present paper in relation to HST imaging, should be car- 
ried out for essentially all imaging of severely crowded 
fields. 
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Fig. 10. Luminosity function of a simulated stellar field with 
n 1 /, ~ 15.8, as observed for the central field, but eener- 
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ated from only one type of stars, namely stars of magnitude 
rhi = 22.78 (vertical line), see text. Note the similarities to the 
observed luminosity function (Fig. 7) . The median indicated at 
the top is 21™ 5. (The histogram above has been normalised to 
the same number of objects as that of Fig. 7) . 
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